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Premixed gas flames in mixtures of CH 4 , 0 2 , N 2 , and C0 2 were studied numerically using detailed 
chemical and radiative emission-absorption models to establish the conditions for which radiatively induced 
extinction limits may exist independent of the system dimensions. It was found that reabsorption of emitted 
radiation led to substantially higher burning velocities and wider extinction limits than calculations using 
opticallv thin radiation models, particularly when C0 2 , a strong absorber, is present in the unburned gas. 
TVvo heat loss mechanisms that lead to flammability limits even with reabsorption were identified. One is 
ibafr for dry hydrocarbon-air mixtures, because of the differences in the absorption spectra of H 2 0 and 
C0 2 , most of the radiation from product H 2 0 that is emitted in the upstream direction cannot be absorbed 
by the reactants. The second is that the emission spectrum of C0 2 is broader at flame temperatures than 
, ambient temperature; thus, some radiation emitted near the flame front cannot be absorbed by the reac- 
tants even when they are seeded with C0 2 . Via both mechanisms, some net upstream heat loss due to 
radiation will always occur, leading to extinction of sufficiently weak mixtures. Downstream loss has prac- 
tically no influence. Comparison with experiment demonstrates the importance of reabsorption in C0 2 - 
* diluted mixtures. It is concluded that fundamental flammability limits can exist due to radiative heat loss, 
but these limits are strongly dependent on the emission-absorption spectra of the reactant and product 
'gases and their temperature dependence and cannot be predicted using gray-gas or optically thin model 
parameters. Applications to practical flames at high pressure, in large combustion chambers, and with 
exhaust-gas or flue-gas recirculation are discussed. 


Introduction 


nmability limit studies are important for as- 
j of fire safety in many environments and for 
ermining the operation limits of combustion de- 
3 ' Despite many years of study, the mechanisms 
f flammability limits of premixed gas flames are still 
| well understood. In particular, it has not been 
Mished whether “fundamental” limits exist in- 
ndent of apparatus size and geometry. Many 
\ mechanisms cause apparatus-dependent extinc- 
limits including flame stretch [ 1 , 2 ], buoyant 
section [3,4] and heat losses to cold walls [5-7]. 
HP*. n these losses are eliminated by employing large 
bustion vessels at microgravity, heat loss via gas- 
! radiation is probably the dominant extinction 
tism [ 8 — 11 ]. 

Nevertheless, it is unclear whether radiative loss 
5 fundamental limits because emitted radiation 
1 be reabsorbed, a factor not considered by the 
ones [6,7] and computations [10,11]. Whether 
jpebsorption is important is usually assessed by com- 
Qg the system dimension to the burned-gas 


Planck mean absorption length (L P ), defined as the 
mean monochromatic absorption coefficient (re) 
weighted by the Planck function: 


■t. 


f 


K(u)G(u)du\ 


G(u) = “ — - — 7 ; u * hcca/kT (1) 


7T* e u - 1 


where h , c, cd, k, and T denote Planck’s constant, 
light speed, wavenumber, Boltzman’s constant, and 
temperature, respectively. In the experiments in 
Refs. [8,9], the optical thickness x = X/L ?1 where X 
is an apparatus dimension, was generally small (with 
exceptions noted later), consequently optically thin 
conditions (no reabsorption) applied. Theoretical 
studies of premixed-gas flames seeded with inert 
particles [ 12 ] predict that with reabsorption, net heat 
losses decrease, burning velocities (S L ) increase, and 
the flammability limit equivalence ratio ( 0 O ) de- 
creases compared with values without reabsoiption. 
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Fig. 1. Example of model predictions for <t> = 0.70, 
y « 0.30 (near the lean Ham inability limit for optically thin 
flames) for adiabatic, optically thin, and reabsorbing 
flames. The corresponding burning velocities are 7.9, 5.9, 
and 12.7 cm/s, respectively. The optical thickness t =* L\f 
Lp is 2.39. Upper plot: temperature (T) and radiative power 
(Q R ), —0.5 < x < 1.0 cm; lower plot: temperature only, 
— 20 < x < 20 cm. 


Corresponding microgravitv experiments in particle- 
laden methane-air mixtures [13] support these pre- 
dictions, Gaseous flame behavior should be different 
from that of particle-laden flames because gases emit 
and absorb in spectral bands, whereas particles ex- 
hibit approximately gray-gas behavior. No compu- 
tational studies of gaseous flames with detailed ra- 
diative emission/absorption models have been 
performed to test their effects on S L and 0 O - Con- 
sequendy, our goal is to model premixed-gas flame 
propagation with detailed radiative emission/absorp- 
tion effects and compare results with experiments 
and theoretical predictions. 

In addition to microgravity studies, modeling of 
premixed-gas flames with reabsorption is relevant to 
combustion at high pressures ana in large furnaces 
because frequently r > 1. For example, at 40 atm, 
a typical pressure for premixed-charge internal com- 
bustion engines, r *■ 1 for cylinders of radius 4 cm. 
Similarly, t > 1 in atmospheric-pressure furnaces 
larger than 1.6 m. Moreover, significant amounts of 
absorbing C0 2 and H 2 0 are present in the unbur- 
ned mixtures of combustion devices employing ex- 
haust-gas or flue-gas recirculation. 


Numerical Model 

The energy and chemical species conservation 
equations for steady planar premixed-gas flames 
were solved using a CHEMKIN-based code [14] 
with arc-length continuation [15,16]. This code was 


employed previously to model optically thin radiat- 
ing counterfiow flames [16]. For this studv, radiative 
transport including both emission and absorption 
was computed using the statistical narrowband 
model with exponential-tailed inverse line strength 
distribution [17]. The radiative transfer equations 
were solved for wavenumbers between 150 and 9300 
cm" 1 with 25 cm" 1 resolution using the S6 discrete 
ordinate method. Radiation parameters for CO«, 
H 2 0, and CO were taken from Ref [18]. CH 4 ra- 
diation was not included because the necessarv spec- 
tral data were not available, but CH 4 radiation is 
minimal even for optically thin conditions [11]. Fur- 
ther details of the model, including accuracy consid- 
erations. are presented elsewhere [19]. Methane ox- 
idation was modeled using an 18-species, 58-step 
chemical mechanism [14]. The spatial position (x) = 
0 was defined as the location where T = 325 K for 
radiation-free flames. Except where noted, the up- 
stream and downstream boundary locations were x 
= - L ( = -30 cm and x = +L 2 = +400 cm. 
respectively. Upstream boundary conditions were 
ambient temperature (300 K), and composition with 
inflow velocity S L and downstream boundary con- 
ditions were zero-gradient. Ambient-temperature 
blackbody walls were assumed at both boundaries. 

Mixtures of CH 4 + (0.210 2 + (0.79 — y)N 2 + 
yC0 2 ) were examined, with y varied to assess reab- 
sorption effects by substituting emitting/absorbing 
C0 2 for radiatively inactive N 2 . Ambient H 2 0 was 
not considered because experiments at standard 
conditions allow at most 3% H.?0 without conden- 


Results 

Figure 1 shows typical computed results. For op- 
tically thin conditions, the volumetric radiative loss 
(Qr) ~ 4<r(T 4 - 7$)/Lp, where tx is the Stefan- 
Boitzmann constant and T 0 the ambient tempera- 
ture, is always positive (upper plot). With reabsorp- 
tion, is negative atx < 0 because some radiation 
emitted at high T (x > 0) is reabsorbed at lower T 
(x < 0). This decreases the net loss and preheats the 
unbumed mixture so that, consistent with the theory 
in Ref. [12], the peak temperature (T*) exceeds ad- 
iabatic (T a< x). With optically thin radiative loss, T de- 
cays downstream to T 0 (lower plot). With reabsorp- 
tion, T still decays downstream (but on a much 
longer scale); thus, some net loss still occurs. All 
flames exhibit conventional convective -diffusive 
zones of thickness 6 3 oc/Sl *“0.1 cm, where a i* 
the temperature-averaged thermal diffusivity. The 
reabsorbing flame additionally exhibits a much 
longer upstream convective -radiative zone (lower 
plot) whose length is comparable to L P (*•19 cm)- 

Fig. 2a shows that reabsorption effects in 
mixtures (y = 0) are minor because <£ 0 is reduced 


\ 
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tC. 2. Predicted burning velocities and peak flame tem- 
tures under adiabatic conditions, with optically thin 
Itive losses and with a radiative model including reab- 
tton effects, in this and subsequent figures, the tum- 
JOint limit is shown, but for clarity, the unstable lower 
Horn branch is omitted, (a) Lean CH 4 -air mixtures, (b) 
i CH^O.SICL + 0.49N 2 + O. 3 OCO 2 ) mixtures (y = 


By 3.8% and, even with reabsorption S Ll is always 
bin than the adiabatic burning velocity (S L a j). How- 
is reduced by 25%, and S L for the limit 
teure decreases similarly (26%). For both optically 
Km and reabsorption models, is close to 

P® theoretical prediction e ~ 1/2 [6,7]. A simple pic- 
bt emerges from these results. Radiation from hot 
M) and C0 2 can escape only if it is emitted up- 
frsam, because downstream radiation will be reab- 
j&fced by other HoO and C0 2 molecules. Thus, 
absorption approximately halves the net heat loss. 
P®ory for optically thin flames [6,7] predicts S [ j hn 
■SpR ; thus, if were halved, Su ini should de - 
29%. Consequently the classical quenching 
p ec h a nis m approximately applies in this case. 
•Flame characteristics are quite different with y = 


Wave number (1/cm) 


Fig. 3. Radiative properties of combustion gases as a 
function of wavenumber (<a). (a) Absorption coefficients 
k((o) for pure CCL and H 2 0 at 1300 K and 300 K. (b) 
Radiative heat flux at cold (upstream) boundary for <P — 
0.50. y - 0.30 (r = 2.39) with reabsorption effects in- 
cluded. Note that the units on the vertical axis are W/m 2 
per m~ l = W/m (i.e., the radiant flux per unit wavenum- 
ber). Also shown is te/*,, where l h is the Planck function for 
a blackbody at T* for this flame (i.e., the maximum possible 
flux.) 


0.30 (Fig. 2b). As with y = 0, for optically thin ra- 
diation, S| and T* are lower than adiabatic values, 
but with C0. 2 substitution, because of reabsorption, 
these quantities are significantly higher than adia- 
batic values. The effect on <P {) is substantial (0.682 
for optically thin conditions vs. 0.442 with reabsorp- 
tion). Nevertheless, flammability limits still exist. 

Figure 3 elucidates two mechanisms that cause net 
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Fig. 4. Effect of upstream domain length (L\) on flam- 
mability limit properties: equivalence ratio (<£ 0 ) and burn- 
ing velocity (S Ui J for y = 0.20. The corresponding flam- 
mability limit for optically thin conditions is <Pq = 0.62. 


heat loss and thus flammability limits even with reab- 
sorption, designated I and II. Mechanism I is that 
*(<y) is very different for H. 2 0 and CO* (Fig. 3a); 
because H 2 0 is absent from the unbumed gas, most 
H 2 0 radiation emitted upstream cannot be reabsor- 
bed. Mechanism II is that k(oj) for C0 2 is broader 
at T* than T {) (Fig. 3a); thus, some C0 2 radiation 
emitted upstream from temperatures near T+ can- 
not be reabsorbed farther upstream where T T 0 . 
Mechanism II can also occur for H 2 0 but not for 
the dry reactants studied here. Both mechanisms ap- 
ply for arbitrarily large domains. The manifestation 
of these mechanisms is seen in the spectrally re- 
solved radiative flux at the upstream boundary (x = 
— L t ) (Fig. 3b). The spectrum for C0 2 + H 2 0 + 
CO is similar to the C0 2 4- HoO spectra at 1300 K 
(Fig. 3a) less the C0 2 spectrum at 300 K, indicating 
that losses arise mainly from the differences between 
the bumed-gas emission and unbumed-gas absorp- 
tion. H 2 0 emission comprises most of the loss 
(mechanism I), but some occurs from C0 2 emission 
(mechanism II). 

These observations show that fundamental flam- 
mability limits due to radiative losses will exist in 
gaseous flames no matter how large the domain or 
what absorbing gases are present, because some ra- 
diation emitted from the high-temperature region 
cannot be reabsorbed by the unbumed gases (mech- 
anism II). These limits depend on the temperature 
effects on the absorption spectrum and cannot be 
predicted via simple mean absorption coefficients, 
as optically thin limits can. Of course, if sufficient 
quantities of inert particles, soot, or other quasi- 
graybody absorbers are present, complete reabsorp- 
tion could occur. 

Loss at the downstream boundary (x = L 2 ) is 
much less important for several reasons. A disap- 
pearing reactant (e.g., CH 4 ) can produce some loss 


at x = Lo via mechanism I, but CH 4 radiation was 
not considered, and the loss would be much less than 
the upstream H 2 0 and C0 2 loss because CH 4 dis- 
appears near x = 0 just as T rises to values where 
significant radiation could be emitted. Some down- 
stream loss via mechanism II can occur because of 
the downstream temperature gradient, but the gra- 
dient and total decrease in T are much less than the 
corresponding upstream values, leading to much 
lower downstream loss. A third mechanism of radi- 
ative loss occurs at the downstream boundary be- 
cause of the blackbody wall with T = T () , but due to 
reabsorption, its influence is confined to the adjacent 
region of thickness Lp (**18 cm). Thus, T* and S L 
are unaffected, as was verified by changing L 2 from 
400 to 100 cm. 

Losses via mechanisms I or II can occur onlv for 
wavenumbers where k(o)) is nonzero on one side of 
the flame but changes to zero over a length smaller 
than the scale {/c(a>)) " 1 over which reabsorption oc- 
curs. Changes occur on the scale 6 for temperature 
and D/S l for species i, where D t is the diffusion 
coefficient. Because the Lewis numbers aJD i are 
close to unity, S is an appropriate scale for both spe- 
cies and temperature changes. Hence, these criteria 
become 


\dx± 

dx 


dXi 


\d(l/K(co))\ 


d ln(/c(a>)) 
d HXi) 


«|S 

■ ft 


> k(co)S => k(co)S < 1 


- 


dT 


\d(l/K(aj))\ 


k(w)S 


d InMco)) 
d ln(T) 


( 2 ) 


where Xi * s the mole fraction of species i and the fact 
k — Xi has been used. Fig. 3b shows evidence of 
these criteria. The loss due to H 2 0 mimics at(w) ex- 
cept where tc(co) is very large and reabsorption can 
occur within the convective -diffusive zone where 
£h 2 o changes rapidly. For C0 2 , practically no loss 
occurs where k(co) is large at both 300 K and 1300 
K (Fig. 3a), but substantial loss occurs for w on the 
“wings” of these peaks where dtddT is large. 

The convective-radiative zone at x < 0 has a char- 
acteristic thickness L P S; thus, L } can influence 
reabsorbing flames drastically. Fig. 4 shows that vo 
decreases as L | increases because more reabsorption 
(thus lower net Q R ) occurs with larger L\. Because 
S L .1 im ~ Cr 2 * Su> m also decreases. Significant reab- 
sorption effects occur even for L { = 1 cm U 
0.054) because C0 2 has absorption bands with k(&) 
up to 4000 m _l atm' 1 (40 cm' 1 atm - 1 ) (Fig- 
For the limit conditions at L| — I cm, Xu 2 o “ * 

thus, significant absorption occurs on the scaw 
(0.190 X 40 cm' 1 )' 1 * 0.13 cm. This estimated 
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ities under adiabatic conditions, with optically thin ra- 
losses, and with a radiative model including reab- 
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um reabsorption scale is also seen in Fig. 1, 
shows the negative (absorption) portion of 
; on a similar scale. Still, L\ > 100 cm (t > 
is required to obtain domain-independent re- 
because the band "wings ’ have much smaller 
) and thus longer absorption lengths. Conse- 
provides only limited insight on flame 
rties with reabsorption. 
urate determination of Q R for small k is diffi- 
because for each 25 cm ~ 1 spectral window, k(co) 
cterizes the sum of many individual absorption 
of widely varying strengths. For some inaivid- 
but not others, absorption coefficients are 
of 10 30 higher at 1300 K than at 300 K. Some 
overlap in co-space while others are separated 
nt gaps that are required for net loss via 
ms I or II to occur. Clearly, a single k(cd, 
each spectral window cannot capture all rele- 
-“tails. Thus, although the reabsorption effects 
e characteristics proposed here are believed 
quantitative predictions will vary slightly de- 
g on spectral resolution, line-shape model, 
on, especially for large L { . 

5 shows that, for all radiation models, sub- 
C0 2 for N 2 has only minor effects on S L in 
1.0 mixtures, except when y is high enough to 
e I *1 and thus, S L ^ appreciably. In contrast, 
% C0 2 substitution nearly doubles S L for reab- 

ig 0 

= 0.5 mixtures. There are two reasons 
jeabsorption effects are stronger at 0 — 0.5. 

= 0.5 mixtures are much closer to the ra- 
extinction limit (note, for optically thin 0 = 
res, very small amounts of C0 2 cause ex- 
ent due to increased Q R ) and thus benefit 
from the reduction in via reabsorption. 


1 


Second, 0 — 0.5 mixtures have much higher Boltz- 
mann numbers (B), which is a scaled ratio of black- 
bodv emissive power at to total heat release rate 
and thus measures the potential for radiative pre- 
heating to increase S L . For nonscattering media [12] 


_ WTU ~ fl) E 

^ad 


(3) 


where /?, £, R, p ih and Cp represent nondimensional 
activation energy, overall activation energy, gas con- 
stant, ambient density and specific heat, respectively. 
Values of B for y — 0 mixtures are about 11.3 and 
127 for 0—1.0 and 0.50, respectively; thus, reab- 
sorption can increase S L much more in 0 — 0.5 
mixtures. 

Figure 6a shows that, for optically thin conditions, 
C0 2 substitution increases 0 () and S L j im , because 
the additional radiating C0 2 increases Q R . With 
reabsorption, small amounts of C0 2 substitution ac- 
tually decrease 0 O and S L u m due to greatly reduced 
Q r , whereas larger amounts increase 0 O slightly due 
to reduced Fig. 6b shows that for optically thin 
conditions, S uliltT /S Uii a is always close to the theo- 
retical prediction [6,7] e ‘ 1/2 , whereas with reabsorp- 
tion can be greater than 20. 

Figure 7 shows comparisons of our computed ef- 
fects of reabsorption on S L to the theoretical predic- 
tion of Ref. [12] S l /Sl . adl n (Si/$ L *j) = ® f° r 
gases. The theory requires that (T* — TidVTid a 
small, O(Ufi) quantity, which is justified for our con- 
ditions (Fig. 2b). Agreement is poor because the 
computation allows substantial net spectral heat loss 
(Fig. 3b), which is not considered by the theory. 
Agreement is improved with H 2 0 and CO radiation 
suppressed (mechanism I eliminated) and still better 
with temperature broadening of k(co) also sup- 
pressed (mechanism II eliminated) so that the flame 
is practically adiabatic. Still, Si/Sl.*! is much lower 
than theoretical predictions because the theory as- 
sumes graybody radiation, whereas gases radiate 
only in certain spectral bands and thus accelerate S L 
less than graybody radiators would. When gray-gas 
conditions are applied (k(cd) = constant = 26 m~ 1 
= Lp l ), the agreement is more satisfactory, though 
calculated results are now above theoretical predic- 
tions, probably because the theory assumes constant 
thermodynamic and transport properties and single- 
step chemistry. Thus, our model is believed to cap- 
ture the essential elements of flame propagation 
with graybody reabsorption plus the nuances of 
spectrally radiating gases. 


Comparison with Experiment 

For CH 4 -air flames (Fig. 2a), comparison between 
computation and microgravity experiments in spher- 
ically expanding flames [20,21] and tubes [22] and 
earth-gravity experiments specially designed for low 
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S L [23] is more favorable without re absorption be- 
cause optically thin conditions are applied to these 
experiments. The discrepancy between computation 
and experiment is significant even though the chem- 
istry and transport models employed predict S L away 
from the limits quite well [24]. There is other evi- 
dence that different chemical models capable of pre- 
dicting burning velocities of near-stoichiometric 
flames yield varying predictions of near-limit prop- 
erties [25]. 

The only relevant S L data that may exhibit opti- 
cally thick conditions are microgravitv spherically ex- 
panding flame experiments [9,26] in (CH 4 + 40 2 )- 
COo mixtures and counterflow twin-flame 
experiments [27] in CH 4 -(0.21O 2 4- 0.79CO 2 ) 
mixtures (y = 0.79). Neither configuration corre- 
sponds to our planar flames, but comparisons are 
nonetheless considered useful. To compare with 
Refs. [9,26], we chose L\ = L* = 6 cm. which cor- 
responds to a flame radius half the vessel radius. Bv 
this radius, S u was steady, yet the pressure rise in 
the constant-volume vessel was small (<10%). To 
compare with Ref. [27], we chose L y = L 2 = 0.35 
cm, which corresponds to twin flames located mid- 
way between the nozzle exits and stagnation plane 
for 1.4-cm nozzle separation. For both experiments, 
ambient values of L P are about 5.4 cm, but even for 
L\ = Lo ~ 0.35 cm, some reabsorption is antici- 
pated (Fig. 3a). 

For the configuration of Refs. [9,26], the optically 
thin model clearly overpredicts the limit fuel con- 
centration (by 13%) and S uliin (by 350%); thus, reab- 
sorption extends this limit. With reabsorption, Su^ 
is predicted well, indicating that the net loss is cor- 
rectly predicted, but the limit fuel concentration is 
underpredicted (by 16%), perhaps because in spher- 
ical geometry the radiative flux divergence causes 
less radiative preheating than would occur in planar 
geometry. For the configuration of Ref. [27], no nu- 
merical solutions could be obtained for adiabatic or 
optically thin conditions with L { = Lo = 0.35 cm; 
this might explain why in Ref. [27] numerical pre- 
dictions of S L could not be obtained. In Ref. [27], 
results were extrapolated to zero strain to estimate 
the planar S L . Although the accuracy of this method 
for reabsorbing conditions has not been established, 
with reabsorption the computed and experimental 
values of S L agree moderately well (Fig. 8). There- 
fore, even for these small-scale flames, reabsorption 
effects may be important. 


Conclusions 

The effects of radiation on pre mixed-gas flames 
were studied using a detailed emission-absorption 
model for H 2 0, C0 2 , and CO. It was found that 
while reabsorption can increase burning velocities 
and extend flammability Urn its considerably, funda- 
mental limits independent of the system dimensions 



Optical thickness 



Fig. 6. Effect of substitution of C0 2 for N 2 on flam- 
mability limit properties with optically thin radiative losses 
and with a radiative model including reabsorption effects, 
(a) Equivalence ratio (#o). burning velocity (Su^), and 
peak flame temperature (T*). (b) Ratio of to the 
burning velocity of an adiabatic flame in the same mixture 
(SLad). 


exist due to the nature of gas radiation, specifically, 
(1) differences between the spectral characteristics 
of reactants and products and (2) temperature 
broadening of the emission/absorption spectra. The 
results agree well qualitatively and in some cases 
quantitatively with theory and experiments. 

In future work, we will examine stationary ‘flame 
balls,” since modeling of recent space experiments 
[25] suggests dominant reabsorption effects in some 
cases. The spherically expanding flame configuration 
will be studied to compare with the microgravity ex- 
periments already cited [9,26]. The effects of ele- 
vated pressures will be examined because of their 
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-FlG. 7. Comparison of theoretical predictions [12] of 
Jjeabsorption effects on S| to numerical results with all ra- 
jjdtotion included, with CCL radiation onlv, with CO : radi- 
i Btion only with the temperature broadening of the absorp- 
spectrum artificially suppressed, and with Ma>) = 
nt = 26 m ~ 1 to simulate a gray gas. 
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^G. 8. Comparison of numerical predictions of S L with 
eriments on microgravitv spherically expanding flames 
+ 4O 2 )-C02 mixtures [9] and counterflow twin- 
5 experiments in CH 4 -air mixtures with y - 0.79 (all 
i replaced by C0 2 ) [27], 


nee to the internal combuston engine and ga 
bine flames. Collisional broadening of the absorp 
* set rum will likelv be important in these cases 
„ , the effects of exhaust-gas or flue-gas recir 
ttion will be examined via computations at ele 
d temperatures with C0 2 and FLO addition. 

Acknowledgments 

^PDR acknowledges support by NASA- Lewis gran 
*AG3-1523 and NAG3-2124. Dr. Fengshan Liu (Nationa 


Research Council. Canada) contributed to building the 
computational code. 


REFERENCES 

L Clavin. P., Prog. Energy Comb. Sri. 11:1-59 (1985). 

2. Williams, F. A., Combustion Theory , 2nd ed., Benja- 
min-Cummins, Menlo Park. CA. 1985. 

3. Jarosinskv, J., Strehlow, R. A., and Azarbarzin. A., in 
yiineteenth Symposium i Internationale on Combus- 
tion , Combustion Institute. Pittsburgh, 1982, pp. 
1549-1557. 

4. Buckmaster, J. D. and Mikolaitis, D , Coinbust. Flame 
45:109-119 (1982). 

5. Jarosinskv, J., Combust. Flame 50:167-175 (1983). 

6. Buckmaster, J. D., Combust. Flame 26:151-162 
(1976). 

7. Joulin, G.. Clavin, P, Acta Astronautica 3:223-240 
(1976). 

8. Ronney, P. D., in Twenty- Second Symposium (Inter- 
national) on Combustion, Combustion Institute, Pitts- 
burgh, 1988, pp. 1615-1623. 

9. Abbud-Madrid, A. and Ronney, P D., in Twenty-Third 
Symposium (International) on Combustion, Combus- 
tion Institute, Pittsburgh, 1990. pp. 423—431. 

10. Lakshmisha, K. N., Paul P. J.. and Mukunda. H. S., in 
Twenty-Third Symposium { International) on Combus- 
tion, Combustion Institute, Pittsburgh, 1990, pp. 433- 
440. 

1L Law, C. K. and Egolfopoulos, F. N., in Twenty-Fifth 
Symposium ( International ) on Combustion, Combus- 
tion Institute, Pittsburgh, 1994, pp. 137-144. 

12. Joulin, G. and Deshaies, B., Combtist. Sci. Tech. 
47:299-315 (1986). 

13. Abbud-Madrid, A. and Ronnev, P. D.. A/AA /. 
31:2179-2181 (1993). 

14. Kee, R. J., Dixon-Lewis, G„ Wamatz. J., Coltrin, 
M. E., and Miller, J. A., Sandia report SAND86-8246. 

15. Giovangigli, V. and Smooke, M. D., Combust. Sci. 
Tech. 53:23-49 (1987). 

16. Ju, Y., Guo, H., Maruta. K., and Liu, F.,/. Fluid Mech. 
342:315-334 (1997). 

17. Malkmus, W.,/. Opt Soc. Am. 57:3234329 (1967). 

18. Soufinani, A. and Taine, J., Int. ]. Heat Mass Trans. 
40:987-991 (1997). 

19. Liu. F., Gulder, O. L., Smallwood, G. J. ( and Ju, Y., Int. 
J. Heat Mass Trans., 41:2227-2236 (1998). 

20. Ronney, P. D. and Wachman. H. Y., Combust Flame 
62:107-119 (1985). 

21. Strehlow. R. A. and Reuss, D. L„ in Progress in Aer- 
onautics and Astronautics , vol. 73. A1AA. New York, 
1981, pp. 61-89. 

22. Okajima. S., Iinuma. K., and Kumagai, S., in Twentieth 
Symposium ( International ) on Combustion , Combus- 
tion Institute, Pittsburgh, 1984. pp. 1951-1958. 

23. Badami, G. N. and Edgerton, A. E., Proc. Roy. Soc. 

I London) A228:297-422 ( 1955). 








2626 


MICROGRAVITY COMBUSTION 


24. Egolfopoulos, F. N, Cho. P., and Law. C. K„ Combust 
Flame 76:375-091 (1989). 

25. Wu. M. S. and Ronnev. P D. ( in Twenty- Seventh Sym- 
posium ( International ) on Combustion , Combustion 
Institute, Pittsburgh, 1998, pp. 2543-2550. 


26. Abbud-Madrid, A., M. S. thesis. Princeton Univp«iK 
1990. 

27. Zhu. D. L.. Egolfopolous. F. N., and Law. C K j n 

Twenty- Second Symposium ( International ) on Com- 
bustion, Combustion Institute, Pittsburgh 1988 
1537-1545. ’ ' PP 





